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Abstract 

Using a schematic solvable many-body Hamiltonian, one studies a new type 
of proton-neutron excitations within a time dependent variational approach. 
Classical equations of motion are linearized and subsequently solved analyt¬ 
ically. The harmonic state energy is compared with the energy of the first 
excited state provided by diagonalization as well as with the energies obtained 
by a renormalized RPA and a boson expansion procedure. The new collective 
mode describes a wobbling motion, in the space of isospin, and collapses for 
a particle-particle interaction strength which is much larger than the physical 
value. A suggestion for the description of the system in the second nuclear 
phase is made. We identified the transition operators which might excite the 

new mode from the ground state. 
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I. INTRODUCTION 


One of the most exciting subject in theoretical nuclear physics is the double beta decay, 
especially clue to the neutrino-less (Ou/3/3) process [1-4]. Indeed, its discovery would answer 
a fundamental question whether neutrino is a Majorana or a Dirac particle. The theories 
devoted to the description of this process suffer of the lack of reliable tests for the nuclear 
matrix elements. O possibility to overcome such difficulties would be to use the matrix 
elements which describe realistically the rate of 2 u/3/3 decay. In this context many theoretical 
work have been focussed on 2 u/3/3 process. Most formalisms are based on the proton-neutron 
quasiparticle random phase approximation (pnQRPA) which includes the particle-particle 
( pp ) channel in the two body interaction. Since such an interaction is not considered in the 
mean held equations the approach fails at a critical value of the interaction strength, g pp . 
Before this value is reached, the Gamow-Teller transition amplitude ( Mqt ) is decreasing 
rapidly and after a short interval is becoming equal to zero. The experimental data for 
this amplitude is reached for a value of g pp close to that one which vanishes Mqt and 
also close to the critical value. Along the time, the instability of the pnQRPA ground 
state was considered in different approaches. The first formalism devoted to this feature 
includes anharmonicities through the boson expansion technique [5-8]. Another method 
is the renormalized pnQRPA procedure (pnQRRPA) [9] which keeps the harmonic picture 
but the actual boson is renormalized by effects coming from the terms of the commutators 
algebra, which are not taken into account in the standard pnQRPA approach. 

In a previous paper [10] we have proved that the pnQRRPA procedure does not include 
the additional effects in an consistent way. Indeed, if the commutators of two quasiparticle 
operators involves the average of monopole terms then these terms should be considered 
also in the commutators of the scattering terms. If one does so, new degrees of freedom 
are switched on and a new pnQRPA boson can be defined. This contains, besides the 
standard two quasiparticle operators, the proton-neutron quasiparticles scattering terms. If 
the amplitude of the scattering term is dominant comparing it to the other amplitudes, the 
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pnQRRPA phonon describes a new nuclear state. 

The aim of this paper is to show that such a mode appears in a natural way within a time 
dependent treatment. The present approach points out new properties of the new proton- 
neutron collective mode. We use a schematic many-body Hamiltonian which for a single 
j-shell is exactly solvable. In this way the approximations might be judged by comparing 
the predictions of the actual model with the corresponding exact results. Since the semi- 
classical treatment is the proper way to determine the mean field, one expects that the 
present approach is suitable to account for ground state correlations in a consistent way and 
therefore some of the drawbacks mentioned in a previous publication [10], like the breaking 
down of the fully renormalized RPA before the standard RPA breaks down, are removed. 
To understand better the virtues of the present model we compare its predictions with the 
results obtained in a renormalized RPA approach and a boson expansion formalism. Since 
the semi-classical methods have, sometimes, intuitive grounds we aim at obtaining a clear 
interpretation for the new proton-neutron mode, ft is well known that the breaking down 
of the RPA approach is associated to a phase transition. In this respect the semi-classical 
formalism is a suitable framework to define the nuclear phases which are bridged by the 
Goldstone mode. Above arguments justify our option for a semi-classical treatment and 
also sketch a set of expectations. 

This project is achieved according to the following plan: In Section II, we describe the 
model Hamiltonian. The main features of the fully renormalized RPA approach, presented in 
a previous paper, are briefly reviewed. A time dependent variational principle is formulated 
in connection with a truncated quasiparticle Hamiltonian, in Section III. This Hamilto¬ 
nian is the term of the model Hamiltonian which determines the equations of motion for 
the quasiparticle proton-neutron scattering terms, in the de-coupling regime. The classical 
equations of motion and their solutions are presented in Section IV. The new pn collective 
mode is alternatively described through the renormalized RPA approach and boson expan¬ 
sion formalism in Section V. Numerical results are analyzed in Section VI while the final 
conclusions are given in Section VII. 


3 



II. THE MODEL HAMILTONIAN. BRIEF REVIEW OF FRN-RPA 


Since we are not going to describe realistically some experimental data but to stress 
on some specific features of a heterogeneous many nucleon system with proton-neutron 
interaction, we consider a schematic Hamiltonian which is very often [11,12] used to study 
the single and double beta Fermi transitions: 
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4. m (c rjm ) denotes the creation (annihilation) of a r(= p,n) nucleon in a spherical shell 
model state \r\nljm) = |rjm) with r taking the values p for protons and n for neutrons, 
respectively. The time reversed state corresponding to \rjm) is | rjm) = (— y~ m \rj — m) 
For what follows it is useful to introduce the quasiparticle (qp) representation, defined 
by the Bogoliubov-Valatin (BY) transformation: 
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which quasi-diagonalizes the first four terms, i.e in the new representation they are replaced 
by a set of independent quasiparticles of energies: 


E r — \J (e T — A r ) 2 + A£. 


(2.3) 


In the new qp representation, the model Hamiltonian, denoted by H q , describes a set of 
independent quasiparticles, interacting among themselves through a two body interaction 
determined by the images of the \ an d Xi terms through the BV transformation. 

Various many-body approaches have been tested by using not the qp image of H but 
another Hamiltonian derived from H by ignoring the scattering qp terms: 
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(2.4) 


(jP n ) ^ l ^nim^iini] 
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m 

and restricting the space of single particle states to a single j-state. Thus, the model 
Hamiltonian contains, besides the terms for the qp independent motion, a two body term 
which is quadratic in the two quasiparticle operators A', A: 

A\jpn) = E4t™a~, A UP n ) = (■ A \jpn)) ] • (2.5) 

m J 

In a previous publication [10], we showed that going beyond the quasiparticle random 
phase approximation (pnQRPA) through a renormalization procedure, a new degree of free¬ 
dom is switched on, which results in having a renormalized pnQRPA boson operator as 
a superposition of the operators A^(jpn), A(jpn) and scattering terms B^(jpn), B(jpn). 
This picture differs from the standard pnQRRPA approach, where the boson operators in¬ 
volve only the operators A t and A, and is conventionally called as fully renormalized RPA 
(frn — RPA). Obviously, when the amplitudes of scattering terms are dominant, one deals 
with a new kind of collective pn excitation. 

In order to define clearly the distinct features of the new proton-neutron (pn) mode 
revealed in the present paper a brief description of the results obtained in a previous publi¬ 
cation [10] is necessary. The equations of motion associated to the many-body Hamiltonian, 
written in terms of quasiparticlc operators, are determined by the commutators algebra of 
the two quasiparticlc (Af A) and scattering (B\ B) operators defined by eqs. (2.5) and (2.4) 
respectively. Within the f rn — RPA , the exact commutators are approximated as follows: 

A (jpn),A\jpn)\ = Cjpl, 

B (jpn ), B j ( jpn ) ] = Cf p i, 

A (jpn), B\jpn)\ = [A(jpn), B(jpn)\ = 0. (2.6) 

The terms cj^, Cf A n appearing in the r.h.s. of the above equations are the averages of 
the corresponding exact commutators, on the correlated ground state |0 >: 
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(2.7) 
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with Nj T standing for the r (=p,n) quasiparticle number operator in the shell j. The nor¬ 
malized operators 
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satisfy bosonic commutation relations and thereby their equations of motion are linear: 
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The matrix T 3J> depends on the U and V coefficients as well as on the strengths y, xi of the 
two body interactions. The frn — RPA approach dehnes a linear combination of the basic 
operators A^jpn), A(jpn), B^(jpn), A(jpn), 


r f = X] [ x U) A Hj'P n ) + Z{j)D ] {jpn) - Y(j)A(jpn) - W(j)D(jpn )] , 


( 2 . 10 ) 


3 

so that the following commutation relations with its hermitian conjugate operator and the 
model Hamiltonian hold: 
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( 2 . 12 ) 


The operators D\jpn ) are identical with B\jpn) or B(jpn) depending on whether the sign 
of Cjpl is plus or minus. The equation (2.12) provides a set of homogeneous equations- 
called the frn — RPA equations- for the amplitudes A", Y. Z. Vb: 
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while the equation (2.11) yields the normalization equation 


£(* 2 (j) + z 2 U) - Y 2 (j) - W 2 (j)) = 1. (2.14) 

3 

The frn — RPA matrices depend on the renormalization constants C ^which, at their 
turn, depend on the phonon amplitudes. Therefore, the equations (2.12) and (2.7) should 
be self-consistently solved. 

In ref. [10] the frn — RPA equations have been solved both for a proton-neutron dipole- 
dipole interaction, needed for the description of the double beta Gamow-Teller decay and for 
a proton-neutron monopole-monopole interaction used in the calculation of the rates of the 
double beta Fermi decay. Equations obtained in the two cases have some common features 
which, for what follows, are worth being enumerated. 

1) The dimension of the f rn — RPA matrix is twice as large as that of the standard RPA 
and consequently new solutions show up. 

2) The solutions characterized by that the largest phonon amplitude is of type Z define a 
new class of proton-neutron excitations. 

3) Due to the attractive character of the two body interaction in the particle-particle (pp) 
channel, the lowest new state has an energy which is smaller than the minimal absolute 
value of the relative energy of the proton and neutron quasiparticle partner states, related 
by the operators (jpn), B(jpn). 

4) For the N=Z nuclei, this minimal value is vanishing and therefore the lowest mode becomes 
spurious or in other words saying a new symmetry is open. The new symmetry corresponds 
to the restriction Cff = 0, i. e. the average of the third component of the isospin operator 
is vanishing. This means that the system is invariant to rotations around any axes in the 
(. X , Y) plane of the isospin space associated to the (jpn) orbits. 

5) Important quantitative effects are expected for heavy nuclei having the proton and neutron 
Fermi energies lying far apart from each other. 

6) The presence of the additional states influences also the structure of the states lying close 
to those predicted by the standard RPA. Indeed, the actual normalization condition for the 
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phonon amplitudes implies new values for the X and Y weights. Consequently, the strengths 
for /3~ and j3 + transitions are shared by the ”old”-lying close to the standard RPA states- 
and the “new” states- for which the amplitudes Z are dominant. 

7) The standard RPA approach is based on the quasi-boson approximation and therefore 
it ignores some important dynamic effects (only the terms A'A', A'A, AA are considered in 
an approximative manner) and moreover the Pauli principle is violated. By contrast, within 
the frn — RPA all the terms of the model Hamiltonian are taken into account. Also the 
Pauli principle is, to a certain extent, restored. Due to this feature, large corrections to the 
double beta transition amplitude as well as to the Ikeda sum rule are expected by changing 
the RPA to the frn-RPA. 

8) The equations of motion for the A\A and B\ B operators are coupled by the terms 
A'B\ A^B, AB\ AB involved in the quasiparticle Hamiltonian. These terms are multiplied 
by the factors U p V n U p /U n ', U p V n V p ’V n i, V p U n U p >U n > , V p U n V p 'V n ' in the ph — ph interaction (the 
X term) and by U p U n U p 'V n ', U p U n y p 'U n ', V p V n U p 'V n i, V p V n V p 'U n > in the pp — hh interaction 
(the Xi term), respectively. Note that the coupling terms change the number of either proton 
or neutron quasiparticles by two units. The terms bringing the main contribution to the 
equations of motion for the operators A\ A commute with Nj p — Nj n but not with N JP + N ]n . 
By contrary the terms having the dominant contribution to the equations of motion for the 
operators B\ B commute with N ]p + N ]n and not with Nj p — Nj n . None of the two operators, 
Nj p — Nj n , Nj p + Nj n , commutes with the coupling terms. Retaining from the ^-interaction 
the pp — hh terms (those multiplied by U p U n V p ’V n ' ) and from the xi interaction only the 
ph — ph terms (those proportional to U p V n U p >V n ' ) the equations of motion for the operators 
B^, B are decoupled from those for kR and A. One may conclude that the new mode is 
determined by a combined effect coming from the pp — hh and ph — ph terms belonging to 
the y and Xi interactions, respectively. 

9) In the particle representation the frn — RPA phonon operator is a linear superposition 
of ph, hp,pp and hh operators. 

10) In the limit of large pp and negligible ph interactions, the amplitudes Z can be analyti- 



cally calculated. The result is that Z is proportional to either U p V n or V p U n , depending on 
whether the sign of E p — E n is plus or minus, respectively. When this amplitude prevails 
over the other ones, the corresponding mode describes a neutron-hole proton-particle (or a 
proton-hole neutron-particle) excitation of the mother nucleus ( N,Z). Therefore the state 
is associated to the (N — 1,Z + 1) (or to the (N + 1,Z — 1)) nucleus. In this case the 
state might be reached by exciting the ground state through the transition operator c^c n (or 
c n c p), which is typical for the /3~ ( or /3 + ) decay. Since the double beta decay is conceived 
as taking place through two successive (3~ transitions, one expects that this process is also 
influenced by considering this new state as an intermediate state characterizing the odd-odd 
neighboring nucleus. 

11) When the pp interaction is small the amplitude Z is proportional to U p U n if E p > E n or 
to V p V n in the case E p < E n . The new mode characterizes the nucleus ( N + 1, Z + 1) in the 
first case and the nucleus (N — 1, Z — 1) in the second situation. The transition operators 
which could excite these states are obviously of the types cj,<4 and c n c p , respectively. 

Note that the restriction of the phonon operator to the scattering terms resembles the 
standard RPA boson operator written in the particle representation. This comparison has, 
however, only a formal value since in the quasiparticle representation there is no Fermi energy 
and therefore one cannot speak about quasiparticlc-quasihole excitations. Similar features 
are met in solid state physics for the description of electron excitations in narrow energy 
bands, spin waves and plasma oscillations [13]. In nuclear physics, the scattering terms 
have been also considered but not for proton-neutron excitations. Indeed, using the thermal 
response theory, Tanabe [14] studied the charge conserving phonons in nuclear systems at 
a finite temperature. It seems that the contribution of the scattering terms to the charge 
conserving bosons, does not survive at vanishing temperature [15]. Moreover, the dispersion 
relation for the mode energy cannot be obtained from a linearized set of equations as it is 
required by the spirit of the RPA approach. 
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III. SEMI-CLASSICAL TREATMENT 


As we already mentioned, the scope of the present paper is to study the pn mode caused 
by the quasiparticle scattering terms within a semi-classical approach. In this formalism 
the renormalization condition (2.6) is missing and therefore the harmonic motion of the new 
degrees of freedom hinges on a more physical ground. Moreover, we address the question 
whether this mode survives when the non-scattering terms are switched off. Thus, it is 
worth to know if such a mode appears only when the scattering terms accompany the two 
quasipatricle operators or it might be determined by the scattering terms alone. 

From the brief presentation of the previous Section it is clear that the mode does not ap¬ 
pear within the RPA approach. Indeed, it occurred within the frn — RPA after a consistent 
renormalization was performed (i. e. not only the operators A', A where renormalized but 
also B t and B). If that mode is a signature of the higher RPA formalisms, then it should 
also appear within the semi-classical formalism as well as in the boson expansion framework. 
As we shall see the semi-classical approach is able to predict the mode even in the harmonic 
approximation, the mode being associated with the small oscillations of the system around 
a static correlated ground state. Moreover, the semi-classical frame is expected to allow us 
an intuitive interpretation of this new type of excitation. 

We recall that the higher order corrections to the standard RPA approach are frequently 
studied, with different purposes, using a single j case and ignoring the scattering terms. The 
procedure has the advantage that the resulting Hamiltonian is exactly solvable. Therefore 
the quality of the adopted approximations may be tested by comparing the predictions with 
the corresponding exact results. 

To touch the goal of the present paper we adopt a similar point of view. Indeed, if the 
coupling terms (mentioned at the point 8 of the previous section) are ignored, the equations 
of motion for the scattering operators are decoupled. Moreover the motion of these operators 
is determined also by an exactly solvable Hamiltonian, which reads: 

Hj£ = E' p N p + E' n N n + A + A 2 (B t2 (pn) + B 2 (pn)), (3.1) 
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where the following notations have been used: 


E’ r = E„ + ( Xl - x )Vl 

E' n = E n + ( X + Xi)V?(V:;-Ul), 

A, = x(Kul + v 2 v 2 ) - Xl(UX + K 2u n), 

A2 — — (x. + Xi)UpU n V p V n} (3.2) 

N r ^ v Q j Tim® J Tjm' (3.3) 

m 

Also, to simplify the notation we omitted the quantum number j for the operators 
D^(jpn), B(jpn) as well as for the U,V coefficients and quasiparticle energies. This model 
Hamiltonian will be studied within a time dependent variational formalism. Therefore, some 
static and dynamic properties will be described by solving the equations provided by the 
time dependent variational principle (TDVP) 1 : 

= (3.4) 

If the variational state T) spans the whole Hilbert space describing the many-body system, 
solving the equation (3.4) is equivalent to solving the time dependent Schroedinger equation, 
which would be a very difficult task. In the present paper, the trial function is taken as: 

|T) = exp[zB^ (pri) — z* B(pn)]\NT — T ), (3.5) 

where lA^TX^) denotes the common eigenfunction of the quasiparticle total number (N), the 
quasiparticlc isospin squared (T 2 ), and its z-axis projection (T z ) operators, respectively, z 
is a complex function of time and z*, the corresponding complex conjugate function. We 
justify this choice by the symmetry properties of the model Hamiltonian. Indeed, let us 
note first that commutes with the quasiparticle total number operator. Moreover, it 
can be written in terms of the quasiparticlc total number operator and generators of the 
SU(2) isospin algebra 


1 Throughout this paper the units of ft = 1 are used 
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(3.6) 


7+1 = ~-^ B \P n ), 

T-i = ~^B{pn), 

To = ^(Np- Nn). 

Due to this property of H^, the function |'h), which is a coherent state for the SU(2) group, 
is the most suitable for a semi-classical treatment. 

Before closing this section we would like to write the trial function in a form which 
suits better the further purposes. Using the Cambel Hausdorff factorization [16] for the 
exponential function, as explained in Appendix A, one obtains: 

|tf) =Ue aBHpn) \NT -T), 

A/” = (1 + a*a)~ T . (3.7) 

where a depends on the polar coordinates, z = pe l<p : 

a = ta n(p)e*A (3.8) 


IV. EQUATIONS OF MOTION 


In order to write the equations of motion provided by the TDVP (3.4), we need the 
matrix element of as well as of the time derivative operator, These can be evaluated 
by direct calculation, using the expressions (3.5) when the average of is considered and 
(3.7) for the classical action. The result is: 


= -T(E' V - E' n + 2Aj) + j(E' p + E' n ) + 2 T(E'„ - E' n + 2^)-^- 


+2T(2T — 1) A, 


• • * 

, T , d . T . _,cr* a - a a 

(T — T) = T -. 

dt 1 + a* a 


a a 


+ X 


a* 2 + a 2 


(1 + a*a ) 2 (1 + a*a ) 2 


(4.1) 


Considering a, a* as classical phase space coordinates, the TDVP equation (3.4) yields the 
following classical equations of motion, describing the nuclear system: 
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(4.2) 


&H _ . T a 

da ^\l + a*a ) 2 ’ 

dH _ T a 

da* (1 + a*a ) 2 

Here 7i denotes the classical energy function: 




(4.3) 


In order to quantize the classical trajectories satisfying the equations (4.2) as well as to have 
an one to one correspondence between the classical and quanta! behaviors of the nucleon 
system, it is convenient to chose those conjugate variables which bring the equations of 
motion in a canonical Hamilton form. A possible choice of the coordinates with the above 
mentioned property is 

2 T 

r = -. 

1 + a*a 

1 

0 = — —(In cv —In a*) =—tp. (4.5) 

Zj L 



Indeed, in the new variables the classical equations read: 


dH 

aT = -*’ 

m _ . 

50 


(4.6) 


with the classical energy: 


N 

H = T(E' p -E' n + 2\ 1 ) + -(E' l ,+ E ') 

2T — 1 

~( E p~ E 'n +2Ai)r+ 2T r(2r-r)(Ai + 2A 2 cos2i^). 


(4.7) 


Note that r has the significance of a generalized coordinate while 0 that of generalized 
linear momentum. Due to the generalized momentum 0, the equations motion are not 
linear and therefore analytical solutions are not obtainable. The equations can however 
be approximatively solved if they are linearized around the minimum point of the energy 
function: 
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(4.8) 


(2T — l)(Ai — 2A 2 )_ 

In order that the minimum exits, it is necessary that the generalized coordinates satisfy a 
consistency condition, required by the definition range of r: 

0 <r< 2T. (4.9) 

By means of (4.8), this provides a constraint for the strengths of the two body interactions. 
The linearized equations, written in terms of the deviations 

q — r — r, p = ifb— % ;, (4-10) 


are of harmonic type: 


. 2T - 1 

P= — 2— 7^ —(Ai — 2A 2 )g, 

‘l- iNzl f (2T- r)X 2 p. 


(4.11) 


These describe a harmonic motion for the conjugate coordinates, with the angular fre¬ 
quency: 

9T - 1 - „ , 

(4.12) 


w = 2^^[-A 2 (A 1 - 2A 2 ) r (2 T- r)]i 


T 

The condition that a; is a real quantity brings an additional constraint for the strength 
parameters XiXv 


X > Xi 


U p Vn - VpU n 


(4.13) 


\ UpUn + VpVn / 

As we said already before, the schematic model has the advantage, over the realistic for¬ 
malisms, that allows us to compare the approximative solutions with the exact one. For the 
particular Hamiltonian used in the present paper, the exact eigenvalues can be obtained by 
diagonalization in the basis NTM). Indeed, in this basis the model Hamiltonian has the 
following non-vanishing matrix elements: 

1 


{NTM\H$\NTM) = -{E' p + E' n )N + (E' p - E' n )M + A,(T + M)(T - M + 1), 
(NTM + 2\H$\NTM) = A 2 [(T - M - 1 )(T - M)(T + M + 1)(T + M + 2)]^ , 


{NTM\H$\NTM + 2) = (NTM + 2\H$\NTM). 


pn 


(4.14) 
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V. THE RENORMALIZED RPA AND BOSON EXPANSION 


Within the RPA approach, the renormalization of the quasiparticle mean field due to 
the two quasiparticle interactions is usually ignored. Therefore the Hamiltonian considered 
is: 

H qp = E p N p + E n N n + A 1 B\pn)B(pn) + A 2 (H t2 (pn) + B 2 (pn)), (5.1) 

The operators B\ B satisfy the commutation relation: 

[B{pn),B ] {pn)\ = N n - N p . (5.2) 

If the r.h. side of the above equation is replaced by its average on the ground state, 

C=(0\N n -N p \0) (5.3) 


which is to be determined, then the operators B, B^ become bosons, after the following 
renormalization 


B \ pn ) = -^=£ f (pn), B(pn) = ~^=B(pn), 
if C is positive, while for negative C the renormalized operators are: 


~ 1 ~ 1 

B\pn ) = B ( pn ), B(pn) = B'(pn). 


#i 


v'lci 


(5.4) 


(5.5) 


Suppose, for the time being, that C > 0. If that is not the case the corresponding calculations 
can be worked out in a similar way. The equations of motion for the renormalized operators 
are: 

H qp , B\pn)\ = (Ep - E n + XiC)B^(pn) + 2A 2 CB(pn), 

H qp , B(pn) = -2A 2 CB ] (pn) - (E p - E n + X 1 C)B(pn). (5.6) 

Since the equations are linear in B^ (pn) and B(pn ), one can define the phonon operator 

T f = XB ] ( pn) - YB(pn ), (5.7) 
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with the amplitudes determined such that the following equations are fulfilled: 


H qp ,r'} = uT\ 


r,r f 


= l. 


The first equation provides the dispersion equation for the mode energy 


LU = 


(E p - E n + AiC) 2 - 4A \C 2 


while the second one the normalization relation for phonon amplitudes: 


A' 2 - Y 2 = 1. 


The renormalized RPA vacuum is defined by 


T1 0 ) = 0. 


The solution of the above equation is: 


|0> = e-s(x) e 2x B ~ I NT, —T). 


(5.8) 


(5.9) 


(5.10) 


(5.11) 


(5.12) 


Then the renormalization constant C can be exactly evaluated: 

C = 2T -2 + ^-. (5.13) 

A- 

Since T > 1, the constant C is always positive. The equations of motion allow us to express 
the amplitude Y in terms of X: 

Y = ^\n-( E p~ E n + AiC)]X, (5.14) 

which together with the normalization condition (5.10) determines fully the amplitudes X 
and Y in terms of C and c o. Inserting the result for X into the equation (5.13), one obtains 
an equation for C as a function of uj. This and eq.(5.9) form a set of two nonlinear equations 
for the unknowns u and C. 
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As we mentioned before, another way to improve the RPA treatment is to use the boson 
expansion concept. Through this procedure, the SU( 2) algebra, with the fermionic genera¬ 
tors t± i,To defined by eq.(3.6), is mapped to a boson SU( 2) algebra, generated by T±i,To- 
Denoting by b + ,b a pair of boson operators, the SU( 2) algebra generators T±i, Ti can be 
constructed as function of 6 + and b. The resulting expressions are conventionally called as 
the boson expansion of the fermionic generators, respectively. There are three distinct boson 
mappings for the fermionic SU(2) algebra found by Holstein-Primakoff [17], Dyson [18] and 
one of the present authors (A. A. R.) [19], respectively. For the present purpose here we use 
the Holstein-Primakoff (HP) expansion: 


D = -^ + (i - 

To = b + b - T. 


(5.15) 


By a direct calculation it can be checked that, by this mapping, to the operator r 2 it 
corresponds a C-number: 

T 2 = T{T + 1). (5.16) 


The fermion Hamiltonian H qp commutes with the quasiparticle total number and the same 
is true for the generators r ±\, Tq. Therefore the image of the quasiparticle total number 
operator through the HP mapping is invariant against any rotation in the isospin space 
and consequently, according to the above equation, is a C-number. Apart from an additive 
constant, the image of H qp through the HP boson expansion is: 

H® = (E p - E n )f 0 - 2A 1 T +1 T- 1 + 2A 2 (T 2 1 + f\). (5.17) 

Making use of eqs. (5.15), the boson mapping of H qp is a infinite series in the bosons b + , 5, 
due to the square root operators. Expanding the square root operators and truncating the 
result at the second order in bosons, the boson Hamiltonian becomes: 


17 



(5.18) 


= (Ep - En + 2A.TJ6+6 + 2A 2 T(6 +2 + b 2 ). 

For a limited range of the interaction strength, this Hamiltonian can be diagonalized through 
a canonical transformation: 


b + = UB + + VB, 
b=UB + VB + , 

1 = U 2 - V 2 . (5.19) 


The restriction that the ” dangerous” terms have a vanishing strength yields the expression 
for the transformation coefficients and the coefficient, u>i, of the diagonal term B + B\ 



1 

71 


Tl + 


\E n — E„ 


\J (E p — E n + 2XiT) 2 — I 6 A 2 T 2 J 


wi = 


(. E p - E n + 2X 1 T) 2 - 16A \T 2 


(5.20) 


Comparing the expressions of (5.20) and u (5.9), one sees that the two energies are 
identical for the limiting case of X = 1, which is met when A 2 = 0 (see eqs. (5.9) and 
(5.14)). 

At this stage it is worthwhile to make the following remarks: a) When the HP boson 
expansion of the model Hamiltonian is truncated at the second order terms in bosons, the 
quasiparticle total number operator is no longer a C number. Therefore the contribution of 
this term should have been considered in a consistent manner. Moreover the truncation is 
justified only for large values of the total isospin T. b) The same inconsistency appears in 
the calculation of the renormalization constant C. Indeed the expression (5.14) is exact and 
therefore includes all contributions coming from the infinite boson series of the correlated 
ground state, given by (5.14). c) Since the boson mapping (5.15) is an unitary transfor¬ 
mation, the exact eigenvalues of H qp are reproduced by diagonalizing the boson expanded 
Hamiltonian . For the second order truncated Hamiltonian, the canonical transforma¬ 
tion breaks down at a critical value of the attractive interaction strength. However the 
diagonalization procedure is able to find the eigenvalues for any strength of the attractive 
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interaction. The resulting energies exhibits a phase transition (the first derivative has a 
jump) at the critical value of the strength. If the second branch of the energy curve could 
also be approximated by an harmonic mode, describing small oscillations of the classical 
system around a stationary state, this is still an open question [12]. d) The HP boson rep¬ 
resentation provides for the harmonic mode the interpretation of an wobbling motion of the 
system around the total isospin, e) The HP boson expansion is justified (in the sense that 
some eigenvalues of the truncated Hamiltonian are close to the corresponding exact ones) 
when the rotation axis in the isospin space is close to the quantization axis (z axis), which 
is usually taken as the axis to which the maximum “moment of inertia” corresponds. If the 
angle between the rotation axis and z-axis is large the harmonic energy may collapse. In 
this case the quantization axis should be chosen as one of the X and Y axes depending on 
the magnitude of the strength of the rj and Ty terms from the quasiparticle Hamiltonian 
. In this case the boson representation suitable for the low order description should be 
of Dyson type [20]. The harmonic approximation for the new representation describes also 
a wobbling motion of a frequency equal to the square root of the product of the inverse of 
the non-maximal moments of inertia normalized to the inverse of the maximal moment of 
inertia. 


VI. NUMERICAL RESULTS 

The formalism described in the previous sections, has been applied to the case j — 

On the proton level, 6 protons are distributed while in the neutron level, 14 neutrons. Alike 
nucleons interact with each other through pairing forces whose strength are G. p = 0.2MeV 
and G n = 0.4 MeV. From the pairing equations it results the following expression for the 
quasiparticle energy: 

E t = Q = (6.1) 

With the data specified above the result for the quasiparticle energy is: 
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E p — 2 MeV, E n = 1 MeV. 


(6.2) 


According to our previous study, the renormalized RPA ground state involves a small number 
of quasiparticles. For example, for a small strength of the particle-particle interaction, the 
quasiparticle total number is about 2 while for large values of the above mentioned strength 
the number may reach the value 4. Due to this behavior of the correlated ground state we 
considered for the isospin carried by the quasiparticles in the ground state, alternatively the 
values 1 and 2. Although these values vary with increasing the particle particle strength we 
kept them constant. 

The numerical analysis refers to the dependence of the energy u of the new nuclear mode, 
on the strengths of the ph and pp monopole interactions, y, Xi- Aiming at showing how good 
is the semi-classical approach for this new type of pn excitation, we calculated also the exact 
eigenvalues of the model Hamiltonian, by diagonalizing the associated matrix (4.14) within 
the basis \NTM). The results are shown in Fig. 1 and Fig. 2. From Fig. 1, one notices 
that the harmonic mode collapses for a critical value of the attractive interaction Xi- This 
critical value is certainly depending on the repulsive interaction strength. The larger is that 
strength the larger the critical value. In Fig. 1, we have also plotted the normalized energy 
for the first excited state. There are intervals for xi where the energy of the harmonic mode 
approximates reasonably well the exact excitation energy. Moreover, for two values of the 
strength parameter, the exact solutions are precisely reproduced. For the case T=2, the 
two energies, exact and to, are the same for xi — 0 and Xi — 0-3 for y = 1 and y = 0.5 
respectively, but the curves are going apart for the first part of interval and then converge 
to an intersection point close to the critical value. 

The peculiar feature of a; as a function of yi, which distinguishes it from the standard 
RPA modes, consists of its non-monotonic behavior with respect to the increase of the 
strength of the pp interaction. The reason is that in the common cases the mean held is 
constant when the two body interaction is varied, while here by changing xi we change also 
the minimum point for energy and therefore another mean held is obtained. It is interesting 
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to notice that although the ph interaction, the y term, is kept constant, the change of the 
mean field is equivalent to an increase of the effective ph interaction until uj reaches the 
maximum value from where the attractive component of the two body interaction prevails. 

In Fig. 2, the energies u and the normalized energy of the first excited state are shown as 
function of y, the strength parameter of the ph interaction. Both energies are monotonically 
increasing with the increase of the interaction strength. In contrast to what happens in the 
case of Xi dependence, here the change of the mean field by changing the energy minimum 
does not change the repulsive character of the y interaction. The agreement between u and 
the exact energy of the first excited state is reasonable good. 

In Fig. 3. the energies characterizing the harmonic mode predicted by the renormalized 
RPA and semi-classical method are plotted as function of Xi- Also, the exact energy of 
the first excited state is presented. Although they have different trends, the semi-classical 
and renormalized RPA energies are not far from each other for y < 0.45. At the critical 
value y=0.57 the energy yielded by the renormalized RPA is going very fast to zero. This 
behavior is specific to the present model where only the scattering terms are considered. 
Indeed, if the phonon operator includes both the two quasiparticle and scattering terms, 
the corresponding mode collapses for larger Xi- In the semi-classical treatment this happens 
only for very large Xi since the static ground state is changed by increasing Xi- The result 
obtained with the truncated HP boson expanded Hamiltonian (see eq. (5.20)) is very close 
to the result shown in Fig. 3 for the renormalized RPA procedure. 

Comparing the results from Fig. la and Fig. 3, we remark on the following features. 
While the renormalized RPA energy collapses at a relatively small value of xi, the mode 
energy predicted by the semi-classical formalism vanishes for a very large yi, far beyond the 
realistic value, which is Xi — X■ This feature is a consequence of changing the static ground 
state with Xi- The energy behavior provided by the semi-classical method is also different 
from that predicted by the standard renormalized pnQRPA (see for example ref. 9) where 
the mode energy is a monotonic function of Xi and goes asymptotically to zero. 

In this context we recall that the frn — RPA breaks down [10] before the standard 
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RPA does, and that happened due to the fact that the lowest frn — RPA energy is that 
associated with the new collective mode. From the present calculations one sees that this 
is not true within the semi-classical approach and therefore including the scattering terms 
in the expression of the phonon operator does not prevent the treatment of the many-body 
system for a realistic value of the pp-interaction strength. 

The vanishing energies for the new mode, shown in Figs, la and 3 suggest that a phase 
transition occurs according to the corresponding formalisms. As we already mentioned 
this is clearly revealed if one diagonalizes the Hamiltonian given be eq. (5.18) [11,12], In 
the renormalized RPA procedure the new phase is determined by a new minimum of the 
classical energy associated to H^. 2 , reflecting the fact that the A 2 term is the dominant one 
for these values of Xi- In the full-line and dotted-line curves of Fig. la, the corresponding 
energies also vanish at certain critical values which result in having again a phase transition. 
This is reflected in the curve obtained by exact calculations, by the fact that the energy is 
minimum for the critical strength. The increasing branch shown by the exact calculations 
(corresponding to the second nuclear phase) might be semi-classically described by changing 
the trial function, involved in the time dependent variational equations, by rotating it (in the 
isospin space) with an angle which corresponds to the orientation of the axis of maximum 
“moment of inertia”. 

It is remarkable that the far intersection points of the curves obtained by semi-classical 
and exact calculations respectively, are lying close to the critical values of the semi-classical 
description. Also, the first intersection point is not far from the critical value of the renor¬ 
malized RPA treatment. In the classical treatment this feature is well known [21], Indeed, in 
the above quoted reference it is shown, for a triaxial rotor cranked on an arbitrarily oriented 
axis, that for certain critical values of the strength parameters, the period of the harmonic 
orbits is equal to the period characterizing the motion on the closed exact orbit. 
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VII. CONCLUSIONS 


The main result of this paper refers to the existence of an harmonic mode determined 
by the scattering quasiparticle terms, which are usually neglected in the standard RPA 
approach. 

The new mode is described within a time dependent variational formalism with an ex¬ 
actly solvable many-body Hamiltonian. The variational state is a coherent state for the 
underlying symmetry group, which is the SU(2) group. A pair of classical canonical con¬ 
jugate coordinates, which bring the equations of motion to the Hamilton form, is found. 
The classical energy has an interesting structure. It is quadratic in coordinate but highly 
non-linear in the conjugate momentum. Therefore one finds first the stationary point which 
minimizes the energy, and then linearizes the equations of motion around the minimum 
point in the classical phase space. The solution for the linearized equations is harmonic 
and its time period determines the energy of the new mode. Despite the fact the classical 
system has an harmonic motion, the mode does not exist in the standard RPA approach. In 
this sense one may say that the present description corresponds to a ’’renormalized RPA”. 
However as we have seen, by comparing the corresponding predictions, the renormalization 
involved in the semi-classical description is completely different from the renormalization 
described in Section V as well as from the boson expansion method. 

It is known the fact that the topological structure of the energy surface depends on the 
strength parameters involved in the model Hamiltonian. Thus, in the parameters space one 
can define several regions, each of them corresponding to a distinct nuclear phase. Having 
this in mind, we studied the behavior of the new mode energy when the strength parameter 
for the pp interaction (yi) is varied. A particular feature for the semi-classical description is 
that the energy is not monotonic decreasing function of \n but it increases in the first part 
of the interval, reflecting that here the ph two body interaction prevails, reaches a maximum 
value, then decreases and finally vanishes. This property is caused by that for each xi a 
new ground state is determined. This aspect is missing in both the renormalized RPA and 
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boson expansion procedures. 

Since the model Hamiltonian resembles the triaxial rotor which was semi-classically stud¬ 
ied by one of the present authors (A. A. R) in refs. [20,21], the interpretation of the new 
mode is imported from there. Thus, the new mode describes a wobbling motion around a 
given total isospin. 

The vanishing energy is a sign for a phase transition. In the first phase the rotation axis, 
in the isospin space lies close to the z-axis, which has the maximum moment of inertia in the 
region of small Xh while for Xi larger than the critical value (where the energy vanishes), the 
rotation axis lies closer to the (X,Y) plane in the isospin space. While the first phase may 
be described by a HP boson expansion formalism, for the second phase the Dyson boson 
representation is the proper one [20]. In the semi-classical approach, the new phase might 
be described by changing the trial function associated to the first phase, through a rotation 
which brings the z-axis to the actual axis of maximal “moment of inertia”. 

The occurrence of the phase transition can be noticed also in the curve showing the 
exact first excitation energy as function of xi- Indeed at the critical value of xi, this curve 
exhibits a minimum. 

Another critical values of xi are those where the mode energy is equal to the exact exci¬ 
tation energy produced by the diagonalization procedure. For these values the linearization 
does not affect at all the period of the exact closed classical orbit. As a matter of fact, for 
Xi lying close to these points the linearization are best justified. It is interesting to notice 
that these values of Xi he however close to the values where the phase transitions in the 
semi-classical treatment (the far intersection point) and the renormalized RPA approach 
(the near intersection point) take place. This observation allows us to conclude that the 
semi-classical approach works very well for the values of Xi where the renormalized RPA 
breaks down and that the interval where the linearization procedure does not work, ending 
with the critical value where the semi-classical energy vanishes, is very narrow. 

The energy of the new mode vanishes for a value of Xi which is far beyond the physical 
value (xi = x)- I 11 this way the drawback of the f rn — RPA, of breaking down earlier than 
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the standard RPA does, is removed. 

How could the new state be populated? We identified the transition operators which 
could excite the new state from the ground state. The conclusion is that these state can be 
seen either in a f3~ (or j3 + ) decay or in a deuteron transfer reaction experiment. 

The coupling of this mode to other collective states will be studied in a subsequent paper 
using a realistic interaction and a large model space for the single particle motion. 
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VIII. APPENDIX A 

Here to derive the factorization of the trial function T). To this purpose we address the 
following more general question. Which are the t-functions A(t ), B(t ), C(t) satisfying the 
equation 

ghb a pm a nm—Z* < l lim a pm )] _ e A(t)ala ne C(t)(Np-N n ) e B(t)a ] n a p (g 

with the initial conditions 

A(0) = B( 0) = C(0) = 0 (8.2) 

and t a real parameter. Once we solve this problem the needed factorization is obtain from 
(4.1) for t — 1. Taking the first derivative of the eq.(6.1), with respect to t, and identify¬ 
ing the coefficients of the similar operators one obtains the following system of differential 
equations for the three unknown functions, A(t ), B (f), C[t) : 

z = A -2 A C ~ B A 2 e~ 2C(t \ 
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0 = c + B Ae~ 2C{t \ 

-z =Be- 2 °U. (8.3) 

Eliminating the functions B, C from these equations, one obtains the following equation for 
A(t). 

z=A -A 2 z* (8.4) 

which admits the solution: 

A(t) = tan (pt)e t<p . (8.5) 

Here the polar coordinates (p, <p)(z = pe ip ) have been used. Inserting the result for A(t) in 

the eq. (6.3), the equations for the remaining functions can be easily integrated. The result 
is: 

C(t) = — In (cos (pt)) 

B(t ) = tan(pt)e _^¥, . (8.6) 

For the sake of simplifying the writing, hereafter the following notation will be used: 

a = H(l) (8.7) 

Using these results the trial function can be written as: 

| T) = e -2C(l)T e A(l)Bt (pn )| ArT _ e A(\)B^(jpn)\ NT _ ^ ( 8 . 8 ) 

where J\f denotes the normalization factor: 

M = e ~ 2C{l)T = e 2l ^ cosp)T = (1 + MT r (8-9) 
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X 1 [MeV] 


FIG. 1. The energy of the harmonic mode given by the eq. (4.12) and the energy of the first 
excited state normalized to the ground state are plotted as function of xi for T = 1 (a)) and T = 2 
(b)) and several values of the particle-hole interaction strength, y. 
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FIG. 2. The energy of the harmonic mode given by the eq. (4.12) and the energy of the first 
excited state normalized to the ground state are plotted as function of xi for T = 1, xi = 2.3MeV 
(a)) and T = 2,xi = 0.0 MeV (b)). 
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FIG. 3. The energies yielded by the renormalized RPA description (5.9)(full line), semi-classical 
approach (4.12) (dotted line) and by diagonalization procedure (dashed line) are plotted as function 
of xi for T = 1 and x = 3-5 MeV. 
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